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Abstract 

We consider a sum rule for heavy-to-light form factors in soft-collinear effective theory (SCET). 
Using the correlation function given by the time-ordered product of a heavy-to-light current and 
its hermitian conjugate, the heavy-to-light soft form factor Qp can be related to the leading-order 
B meson shape function. Using the scaling behavior of the heavy-to-light form factor in AQco/^b) 
we put a constraint on the behavior of the B meson shape function near the endpoint. We employ 
the sum rule to estimate the size of C,p with the model for the shape function and find that it 
ranges from 0.01 to 0.07. 
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Information on heavy-to-light form factors is important to extract Standard Model pa- 
rameters from experimental results. For instance, the B — > 7r form factor is needed to extract 
V u b or the CKM angle 7 from exclusive B decays |l| . Due to the non-perturbative nature of 
the form factor, techniques beyond perturbation theory need to be employed to try to eval- 
uate these functions. Unquenched lattice results are starting to be available for the B — > n 
form factor |2|. Due to the uncertainties in the lattice technique, the pion is restricted to 
have energies E w ^ 1 GeV, just in the region where the experimental uncertainty is largest. 
Another method for determining the form factor is light-cone sum rules [^J. In Ref. ji| (see 
also 0), light-cone sum rules were investigated using soft-collinear effective theory (SCET) 
0], where it is argued that using SCET allows for a more consistent scheme to calculate 
both factorizable and non-factorizable contributions than the traditional light-cone sum rule 
approach. 

In this paper, we investigate a novel sum rule using SCET, relating the B — ► 7r form 
factor at small q 2 to the B meson shape function which describes the motion of the b 
quark inside the B meson. The shape function cannot be calculated, but it can be extracted 
from the data j^f . After relating the two non-perturbative functions using the sum rule, we 
use the scaling in Aqcd/^6 of the B — ► tt form factor to put constraints on shape function 
models near the endpoint region. By choosing a model for a shape function which satisfies 
the constraints, we can give model-dependent values for the B — > n form factor at q 2 = 0. 

Let us consider the correlation function in SCETi defined as 

r%) = i J d 4 ze^ z (B\T[4(z)J (0)]\B), (1) 

where the leading-order heavy-to-light current in SCET can be written as 

J (x) = e i(p - m » v) - x £ n WY^h v (x). (2) 

Here p = n-pn* 1 /2+p^_ is the label momentum of the collinear quark £ n , which has momentum 
p/j. = anc i we w jU se t = 0. We are interested in the correlation function under the 

same kinematic condition as the forward scattering amplitude of inclusive B decays in the 
endpoint region. We denote the residual momentum r as the remainder of the sum of large 
momenta of the heavy quark, the collinear quark, and the incoming momentum q, with 

r M = rribV + g M — p M ~ C(A-), n ■ r = n ■ q + m&, n ■ r = r± = 0. (3) 

n ■ q is chosen to be negative, and A is a typical hadronic scale. Since n ■ r is of order A, 
the quark field in the intermediate state becomes hard-collinear with p 2 ~ 0(m&A). If we 
take a cut of the correlation function 11(g), the final state contains a pion or other collinear 
particles which include a hard-collinear quark and an ultrasoft (usoft) quark. The Feynman 
diagram for the correlation function at leading order is schematically shown in Fig. ^ 

The dispersion relation of the correlation function can be written as 

n[( Pfl + ,?] = ir d {PB + q f ( Jf lip ( B+ f\ 2 . , (4) 

tt Jml^o {pb + q) 2 ~ [Pb + q') 2 - «e 
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FIG. 1: Feynman diagram for the correlation function II in the timelike (spacelike) region. The 
double line is a heavy quark field, the solid line is a hard-collinear quark with the residual mo- 
mentum k(k'). The dashed line is a spectator quark in the B meson, which is an usoft field in 
SCETi. 

where p B is the momentum of the B meson. With A = m B — m b , we can write 

(Pb + qf = {m B v + qf = [(m b + A)v + q} 2 

= n-p(n-r + A) + C(A 2 ), (5) 

and similarly (p B + q') 2 = n ■ p(n ■ r' + A) + 0(A 2 ). The pion mass is regarded as ml ~ A 2 . 
In terms of the residual momentum n ■ r, the dispersion relation is 

U(n-r') = - dn-r ^ / . (6) 

71 Jml/n-p-K U ■ T — U ■ T — %e 

In the lower bound of the integral, m 2 /n ■ p can be neglected compared to A. 

Using the optical theorem, the correlation function in the timelike region can be written 

2Imn(r) = J2 [ ^n e ^(^l J o(^)l^)(^l^o(0)|S) = 2 ImIT 9 + 2 Imn c , (7) 

where the superscripts 'S" and 'C mean the contributions from the single-pion state and 
from the continuum state respectively. The imaginary part of the correlation function from 
the single-pion state is given by 

2ImIF 7 = J ^^-(27r) 4 5(^ + g- P7r )(fi|J Vn)K|Jo|5), (8) 

where the matrix element between the B meson and the energetic pion 7r n is expressed in 
term of the heavy-to-light form factor 

(n n \J \B) = n ■ p Cp(aO- (9) 

Therefore the correlation function for the single pion in the timelike region can be written 

as 

2\mU s = 2ti 5{{p B + q) 2 - ml) (n ■ p) 2 &^) ~ 2tt 5((p B + q) 2 ) (n-p) 2 C P (^. (10) 



as 
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Using Eqs. (jU), ©, and (fTUj). the dispersion relation Eq. © can be rewritten as 

n(n-r') = — t ^ — - + -/ , 11 

—A —n-r—ie ix Jw' a w — n ■ r — ze 

where w = n ■ r and w' s = 4m^/n • p — A is related to the onset of the continuum states, 
which we take as so = 4m^. = n ■ p(w' s + A). 

At the parton level, the correlation function U(n ■ r') can be obtained from SCET. From 
the matching between SCETi and SCETn, the time-ordered product of the heavy-to-light 
current and its hermitian conjugate can be expressed in terms of the jet function and the B 
meson shape function in SCETn- Starting from Eq. the correlation function is 

n( n ■ r') = -m B J J dn ■ ke^'-^^^Jpin ■ k) 

x(B v \hY(^)^h(p)\B v ), (12) 
where the jet function Jp with P = n ■ p is defined as 

(0\TW%(z) ?W(0)|0) = ip(n ■ z)5{z ± ) J ^ e - iMz / 2 J P (n ■ k). (13) 

The jet function can be computed perturbatively and at tree level it is given by 

4V&) = — i—. (14) 
n ■ k + le 

The B meson shape function is defined as 

(B v \hY{^)Y^h{o) \B V ) = J dn- le inmz/2 (B v \hY5(n -l-n- id))Y ] h\B v ) 

= j dn- le in - m - z/2 f(n ■ OTr^^ 



2 j dn-le in - ln z/2 f(n-l), (15) 



where the residual momentum n-l should be less than A since the shape function has support 
for n ■ I < A. Because n-k — n-l + n-r' from the phase space integral in Eq. (|12|) . the 
correlation function can be written as 

/■A 

II(n ■ r') = —nip / dn ■ lf(n ■ l)Jp{n ■ I + n ■ r + ie). (16) 

J — oo 

At tree level, the correlation function becomes 

n(n • r') = m B [°° dw f ^~ W = n ' l \ (17) 
J—A w — n • r — ie 

where we have changed the integral variable to w — —n ■ I to allow for an easier comparison 
with Eq. (JIIJ). 
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From Eqs. and (fTTj) . we obtain the tree-level relation 



dw- 



n ■ p( 



— + - 

le it 



dw- 



ImU(w) 



-A w — n ■ r" — le —A — n ■ r' 
Taking a Borel transformation on each side of Eq. (|18|). we obtain 

1 



w — n ■ r 



le 



m B 



dwf(—w)e 



-W/W M — TfT 



n ■ p( 2 P e A/wM + 



7T Jw 



dwlmll(w)e 



-w/w M 



(18) 



(19) 



From this relation, we can relate the heavy-to-light form factor to the B meson shape 
function as 



1 






m B 


n ■ p 




m B 


r r°° 


n ■ p 


J-A 


m B 


r w ' s 


n ■ p 


J-A 



1 



dwf(-w)e- (w+A)/wM - - / dwlmU(w)e~ {w+A)/wM 

A 7T Jw' s 

roc 

dwf(-w)e- {w+A)/wM - / dwf(-w)e~ {w+A)/wM 

Jw' s 

dwf{-w)e~ {w+A)/wM . 
By changing the variable w + A — > w, we finally obtain our main result 



m B 



dwf(A-w) e ~ w/wM , 



(20) 



(21) 



n ■ p Jo 

where the Borel parameter wm is given by M 2 = wmu ■ p ~ m^A with wm ~ Aqcd- The 
other parameter w s is given by w s = w' s + A with s = 4m^ = w s n ■ p. Note that w s is a 
measure of the distance to the continuum and is roughly of order AQ CD /mfe. 

When we expand the integrand in a Taylor series near w = 0, Eq. (|21j) becomes 



m B 
n ■ p Jo 

m B 



dw { /(A) - w 



/'(A) + 



/(A) 



Since /(A) = 0, Eq. 



/(AH - -y- 



/'(A) + 



w M 
/(A) 



+ 



+ 



/"(A) + 



2f(A) , /(A) 



+ 



+ ■ 



r(A) + 2f(A) + /(A) 



A7 



.(22) 



simplifies to 
n-p\ 2 iU 6 



/"(A) + 



2/'(A) 



+ 



(23) 



^From various arguments [ill U| the soft form factor (p scales as (Aqcd /m-b) 3 ^ 2 , which puts 
a constraint on the behavior of the B meson shape function near A, though the functional 
form is unknown. Since the left-hand side scales as (Aqcd/^6) 3 , we can constrain the 
scaling of the terms on the right-hand side. With w s ~ A 2 /rrib, the shape function /(A) and 
its derivatives of /(A) should scale as /(A) ~ A, /'(A) ~ 1/A, /"(A) ~ 1/A 3 , • • •, based 
on Eq. ()22j) . In general the shape function has a width of order A and is normalized to 
one, so it has a typical size of 1/A ~ 1/A. There is an additional factor of 1/A for each 
derivative of the shape function, from which we expect the naive scaling like /(A) ~ 1/A, 
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/'(A) ~ 1/A 2 . The constraint from the scaling behavior of the soft form factor does not 
allow this naive scaling. In fact, /(A) is suppressed by A 2 and /'(A) is suppressed by A 
compared to the naive scaling behavior due to this constraint. Practically we put /(A) = 
since it is suppressed by A 2 compared to the naive scaling behavior, and assume that the 
series in Eq. (I23j) converges rapidly such that we choose the terms with the first derivative. 
However this scaling behavior, especially /'(A) ~ 1/A puts a constraint on the endpoint 
behavior for the shape function models. 

We can estimate the size of (p from Eq. (}2T]) given a model of the shape function. We 
adopt the model of Ref . , which is given as 

ife " if (24) 

at leading order in a s . The constraint that the first derivative scales as 1/Aqcd implies that 
b > 2 for this model. The default choice of parameters in Ref. |l2| is A = 0.63 GeV and 
b = 2.93, consistent with this constraint. With this choice of the parameters, and using 
w s = 4m 2 / m b = 0.016 GeV, % = 0.5 GeV, we obtain 



Cp « O.Ol/^. (25) 
y n ■ p 

On the other hand, if we choose w s as 0.05 GeV, which can be considered as a typical scale 
AQ CD /m b with Aq CD ~ 0.5 GeV, we get 



Cp^ 0.045,/^. (26) 
y n ■ p 



The behavior of (p with respect to w s and b is illustrated in Fig. 121 The soft form factor 
with the input parameter b = 2.93 is smaller than 0.1 in the region w s = [0, 0.1] GeV and 
the central value is (p ~ 0.05 However, when we set the input parameter b = 2.0, the soft 




b 



FIG. 2: Dependence of Cp on the input parameter w s and b. In the first plot, the solid line denotes 
the soft form factor with b = 2.93 and the dashed line with b = 2.00. In the second plot, the solid 
line represents the variation on b with w s = 0.016 GeV and the dashed line with w s = 0.05 GeV. 
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form factor can be enhanced almost three times compared with the case of b = 2.93. This 
can be easily understood when we consider Eqs. (|23j) and (|24j) . In the case with b = 2, the 
shape function in Eq. (|24j) scales as /'(A) ~ 1/A 2 , which makes (p ~ A overestimating the 
scaling behavior. The shape function model in Eq. (}2"4"|) with b ~ 3 is reliable considering the 
scaling behavior of (p near the tail of the shape function. 

We can estimate the size of the soft form factor using the shape function model extracted 
from B — > X s 7 . The central fitted values on the exponential model of the shape function 
in Eq. (J24j) are 

A = 0.545 Gev, Ai = -0.342 GeV 2 , (27) 
where Ai is defined as (B v \h v (n ■ iD us ) 2 h v \B v ) = — Ai/3 and, using the relation 0] 



A 



Ai 



(28) 



duj UJ 2 f(uj) = = -y, 

we can obtain the value of the input parameter b = 2.605 under the choice of Eq. (J27|) . Using 
the fitted shape function from the data on B — > X s 7, the soft form factor at tree level with 
la error is given by 



in ■ p 
m B 



(p 



0.0175: 



-0.0342 
-0.0167 



071 n+ 0540 
u.uf iu_ a0580 



[W, 



0.016 GeV), 
0.05 GeV). 



(29) 



The behavior of (p using this shape function model is shown in Fig. El In the heavy-to-light 
soft form factor such as £p, it is reasonable to choose w s as Am^/mi, since the invariant mass 
squared of the lowest continuum state starts at p 2 ~ 4m 2 . Therefore, with this choice of 




. 15 




w s [GeV] 

. 02 0.04 0.06 0.08 0.1 b L J 

FIG. 3: Tree-level estimate of £p with the fitted input parameters extracted from the B — > X s ^ 
data. The solid line shows the soft form factor with the central fitted parameters (A = 0.545 GeV, 
Ai = —0.342 GeV 2 ). The two dashed lines represent the range of the soft form factor with la error, 
where the upper line is obtained with the input parameters (A = 0.781 GeV, Ai = —1.13 GeV 2 ), 
and the lower line with A = 0.485 GeV, Ai = -0.13 GeV 2 . 
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FIG. 4: Feynman diagrams for the lowest-order contributions to Ili(g) in SCETj. 



Wsi we find C,p{B — > n 
factorization approach 



at tree level can be small contrary to the prediction of the QCD 
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The numerical analysis of the soft form factor (p is based on the tree-level relation in 
Eq. (|2*T]) . There can be a few sources of theoretical uncertainties. For example, there may 
be uncertainties in choosing the Borel parameters w s and wm- Since the forward scattering 
is computed in SCETi, the virtuality of the intermediate state, or the pion state is of order 
M 2 = ■ p ~ mi,A, and % ~ A. And w s is the scale from which the continuum 

states start. We take as w s = Arn^/n ■ p or ~ A 2 jm^ with A w 0.5 GeV, assuming that 
m 2 ~ A 2 , which is obviously an overestimation. This is in contrast to the choice of the 
Borel parameters in Ref. jij, in which they take w s ~ wm ~ 0.2 GeV. Their tree-level 
evaluation of the soft form factor corresponds to (p ~ 0.32, while our evaluation is about 
0.01 (w s = 0.016 GeV) and 0.05 (w s = 0.05 GeV), which is smaller by an order of magnitude. 
We find that the variation of % between 0.2 GeV and 0.6 GeV gives 1% (w s = 0.016 GeV), 
10% (w s = 0.05 GeV). Therefore we can argue that the large uncertainty arises from the 
choice of w s . Indeed, we obtain Cp ~ 0.30 in our evaluation with w s = 0.2 GeV, % = 0.4 
GeV, close to the result in Ref. |4|. 

Let us briefly comment on the higher-order corrections to the tree-level relation in 
Eq. (|2H). First, we have to include the radiative correction to n(g), which include the 
radiative corrections of Jo- m addition, we also have to include the spectator interactions. 
When we consider a hard collinear gluon exchange with the spectator quark, the contribu- 
tions can be separated into two parts, the soft form factor (p and the hard form factor (j, 
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141. The hard 



based on the existence of the spin symmetry of the heavy-to-light current 
spectator contributions for (p satisfying the spin symmetry can be distinguished simply by 
the presence of the leading order current Jo in the time-ordered product 0. In this case 
we need to include the higher order interaction terms of SCETj collinear Lagrangian in the 
correlation function. 
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The spin symmetry breaks down explicitly if sub leading currents with ipf are included 
in the time-ordered product. Therefore, when the hard spectator contributions for (j are 
considered, we have to define a correlation function with the subleading current J\, which 
is given by 

J^x) = e i( P- mbV > x ljtiptWY^h v (x). (30) 

In this case, the intermediate state consists of hard-collinear particles and one of the corre- 
lation functions is given by 

ni(g 2 ) =i I d'ze^J d 4 x J d 4 y(S|T[4( 2! )J 1 (0) 4 -4J(a;)«X^ ) (2/)]|S), (31) 

where CfJ = £, n {'rf~ L + )W r ^<?«s is the ultrasoft-collinear Lagrangian. This correlation 
function contributes to the hard form factor. After integrating out the hard-collinear degrees 
of freedom, we have four-quark operators with h v and q us . As seen in Fig. 01 from the four 
possible contractions of the hard-collinear particles in the intermediate state we obtain the 
operators in SCETn of the form 

0^1+,%) = (hY) a p(lt -n.id)(Y^h) b .(q u X) b p(lt-n-id)(Y^q us ) a , (32) 

2 (/+, /+) = (hY)ta ±aP 8(lt - n ■ id)(Y^h) b ■ (q us Y) b ^afS(lf - n ■ id){Y^q us ) a , 

3 (/+,/+) = {hY) a ^5{lt-n-id){Yq u X-{q us Y) b ^5{lt-n 

4 (l+ It) = (hY) a ^a ±aP 6(li - n ■ id){Yq us ) a ■ (q us Y) b ^af <5(l+ - n ■ id)(y*h) b , 

where ex"' 3 = i[lj_il±]/^i an d a ' & are c °l°r indices. In addition, we also have to consider 
the possibilities that the hard-collinear loop diagrams in Fig. 0] can be matched onto the 
nonlocal time-ordered products such as Fig. El in SCETn- When calculating the hard- 
collinear loop corrections, we are confronted with non-vanishing infrared divergences which 
result from internal propagators with virtualities of order A 2 . In this case we expect the 
infrared divergences will be reproduced by the possible nonlocal time-ordered products in 
SCET n . 

Therefore, by matching onto SCETn, we obtain 

rii = JdltdltJ2J t (n-p,n-rJtJt)(B\O t (lt^t)\B) 

i=l 

+ J d i xdSY. J ^ J k®{B\i:[0^x) ) 0^y)}\B) + --- ) (33) 

i,k 

where J iy J- are the jet functions and ® denotes convolutions of If, If which can appear in 
the operators Of, 0\. The matrix elements between the B meson in Eq. (|3*3*|) are new non- 
perturbative parameters to be included by defining additional subleading shape functions 
of the B meson. Similarly, there are also hard spectator contributions to C,p from other 
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FIG. 5: Examples on the nonlocal time-ordered products in SCETn- The internal loops consist of 
the collinear particles with offshellness A 2 . The dotted vertices denote the four-particle operators 
obtained from the matching respectively by integrating out the hard-collinear gluons (a,d) and 
quark fields (b,c). 

time-ordered products with the subleading collinear Lagrangian, which introduce more non- 
perturbative parameters. However, we have no information on these subleading B meson 
shape functions, and hence we do not further consider higher-order corrections and other 
theoretical uncertainties such as the renormalization scale dependence in this paper. 

We have derived the tree-level relation between the soft form factor for B — > n and the 
leading B meson shape function using the sum rule approach. This is achieved by considering 
the forward scattering amplitude of the heavy-to-light current near the endpoint. We derived 
a constraint on the behavior of the B meson shape function near the endpoint using the 
scaling behavior of the soft form factor, (p ~ (A/m^) 3 / 2 . Since we work at tree level only, 
there may be a lot of uncertainties, and there may also be a model dependence on the form of 
the shape function. But the analysis suggests that the soft form factor can be about 0.01 to 
0.07, considering various inputs from the models and using the experimental data on the B 
meson shape function. The numerical value of the soft form factor is approximately smaller 
by an order of magnitude than other sum rule predictions. Since we have not included 
higher-order corrections in a s , this evaluation may have a large theoretical uncertainty, but 
our analysis favors a small value of the soft form factor. 
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